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Ripping Apart at the Seams: The Network of Stripped Gas 
Surrounding M86 
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1 INTRODUCTION 



Observations of elliptical galaxies in clusters with XMM-Newton 
and Chandra have shown that the cluster environment has a pro- 
found impact on the evolution of galaxies. There are several mech- 
anisms by which galaxies in clusters interact with both neighbor- 
ing galaxies and the surrounding Intracluster Medium (ICM). Most 
galaxy-galaxy interaction mechanisms usually require low-velocity 
approaches to be efficient, and typically result in the stripping of 
both diffuse interstellar gas and a fraction of their stellar popula- 
tions. Galaxy-ICM interactions commonly result in the stripping of 
the gas associated with the galaxy halo. In the most dramatic ex- 
amples, ram pressure stripping can result in tails of both hot and 
cool gas trailing behind the galaxies for distances of ~ 100 kpc or 
longer. The physical processes that operate within these tails, how- 
ever, are not fully understood. While simulations have shown that 
simple analytic estimates are fairly accurate in predic ting the radius 
of the gas disc that undergoes ram pressure stripping dGunn & Got! 
1 1972b . complex physics is required to account for the morphol- 
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We present a new study of the Virgo Cluster galaxies M86, M84, NGC 4338, and NGC 4438 
using a mosaic of five separate pointings with XMM-Newton. Our observations allow for ro- 
bust measurements of the temperature and metallicity structure of each galaxy along with the 
entire ~ 1° region between these galaxies. When combined with multiwavelength observa- 
tions, the data suggest that all four of these galaxies are undergoing ram pressure stripping by 
the Intracluster Medium (ICM). The manner in which the stripped gas trailing the galaxies in- 
teracts with the ICM, however, is observably distinct. Consistent with previous observations, 
M86 is observed to have a long tail of ~ 1 keV gas trailing to the north-west for distances of 
~ 100- 150 kpc. However, a new site of ~ 0.6 keV thermal emission is observed to span to the 
east of M86 in the direction of the disturbed spiral galaxy NGC 4438. This region is spatially 
coincident with filaments of Ha emission, likely originating in a recent collision between the 
two galaxies. We also resolve the thermodynamic structure of stripped ~ 0.6 keV gas to the 
south of M84, suggesting that this galaxy is undergoing both AGN feedback and ram pressure 
stripping simultaneously. These four sites of stripped X-ray gas demonstrate that the nature 
of ram pressure stripping can vary significantly from site to site. 

netic fields jRuszkowski et alj|2012b are all expected to influence 
the morphology of the stripped gas and its evolution after it is de- 
posited into the ICM. The relative contributions of these processes 
remain uncertain. 
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Arguably the ideal site to observe these interaction mecha- 
nisms is the nearby, bright Virgo Cluster. The most massive galax- 
ies in this cluster have been studied in detail at multiple wave- 
lengths, and X-ray observations have shown convincing evidence 
for profound interactions between the galaxies and the ICM. Long 
tails of cool, stripped gas have been obse rved trailing in the wak e 
of Virgo elliptic als including NG C 4552 dMachacek et alj|2006ah . 
and NGC 4472 jKraft et alj20Tlh . 

Observations of the Virgo Cluster have shown that it is com- 
posed of three separate sub-groups: one centered on M87, the 
Brightest Cluster Galaxy (BCG); one centered on M86, an elliptical 
galaxy roughly one degree (~ 300 kpc) to the west of M87; and the 
final sub-group centered on NGC 447 2 (also known as M49) ap- 
proximately 4.5° to th e south of M87 teinggeli et"afll 19871 . fl993l ; 
iBohringer et al.|[l994|) . Spectroscopic studies of M86 have shown 
that it has an absolute blue-shifted velocity of ~ 250kms~' while 
M87 has a red-shifted velocity of ~ 1300 km s~ ] . Give n an average 
temperature of the ambient Virgo ICM of ~ 2.3 keV dUrban et al.l 
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I20T1L corresponding to a sound speed of ~ 700 km s '), M86 must 
be traversing through the Virgo Cluster supersonically with a Mach 
number of at least M > 2. 

In the vicinity of M86, additional galaxies are located within 
approximately ~ 100 kpc of one another in projection. The two 
most massive galaxies, M86 and M84, are separated by a pro- 
jected distance of only 60 kpc. Although these galaxies appear 
similar in optical imaging studies, further studies utilizing optical 
spectroscopy and X-ray observations show that they are interact- 
ing with the ambient Virgo ICM in different manners. Observa- 
tions of M86 with Chand ra and XMM-Newton dFinoguenov et al.l 
I2004 1 Randall et al.l 120081) show that it is undergoing significant 
ram pressure stripping, with diffuse X-ray emitting gas as well 
as cooler material trailing behind the galaxy for distances of ap- 
proximately ~ 150 kpc in projection to the northwest. There is 
only sparse evidenc e for ram pressure stripp ing in observations of 
M84 with Chandra dFinoguenov et al.ll2008l) . but there is clear ev- 
idence for a powerful A GN outburst at the center of this galaxy 
dFinoguenov et al.l |2008|). M84 has an absolute red-shifted v eloc- 
ity of ~ 1000 km s 1 dSmith et alj|2000l . iTrager et al.l l2000h . ap- 
proximately consistent with the measured velocity for M87. The 
region in between M84 and M86 is also claimed to b e the site of 
bright intracluster light emission dRudick et al.ll2010h . Two other 
galaxies are detected in the X-rays nearby: NGC 4388 and NGC 
4438. These spiral galaxies have absolute redsh ifted velocities 
of 2566 km s"' (NGC 4388. Ikines & Gelleill2008h and 71kms ' 
(NGC 4438, lKennev et aljl995h . respectively. 

In this work we present, for the first time, an investigation 
into the thermodynamic structure spanning the entire region across 
these four galaxies. Utilizing the full information of five separate 
XMM-Newton observations, including new observations of M84, 
we are able to investigate the temperature and metallicity struc- 
ture around these galaxies and investigate how interactions of these 
galaxies with each other and the ICM have transformed their gas 
reservoirs. 

This paper is structured as follows: in Section 2 we discuss the 
data sets, their initial processing, and the production of images and 
spectral maps; in Section 3 we present our results; and in Section 
4 we interpret their physical implications. We assume a nominal 
distance of 17 Mpc to the Virgo Cluster, at which 1 arcsecond cor- 
responds to 82 pc, and 1 arcminute to 4.92 kpc. 



2 DATA PROCESSING AND ANALYSIS METHODS 

A total of five separate XMM-Newton pointings were analyzed, 
the details of which are listed in Table Q] The data products were 
cleaned and processed using the XMM-Newton Science Analysis 
System (SAS) (version 1 1. 0), and in particular th e Extended Source 
Analysis Software (ESAS, Isnowden et alj|2008h package included 
with SAS. We have utilized events from the two MOS and PN 
detectors for this study. Calibration products and event lists were 
first produced using the standard SAS routines. Good time intervals 
(GTI's) for each event list were determined using the MOS-FILTER 
and PN-FILTER routines. Four of the five pointings had little time 
contaminated by flares, but the most recent observation (of M84, 
Observation # 0673310101, PI S. Ehlert) was heavily contaminated 
by flaring, with only ~ 40 ks of the original 130 ks MOS exposure 
deemed suitable for analysis. No further evidence for flaring was 
observed in any of the event lists after GTI filtering. 



Table 1. Summary of the five XMM-Newton observations utilized in this 
study. Exposure times are the net exposure after all cleaning and processing 
as described in Section [2] for each detector. The M84 observation (Obs # 
0673310101, denoted by an asterisk) was heavily contaminated by flares. 
The original event list was split into two separate event lists for each MOS 
instrument, the summed exposure time of which is shown below. 



Obs# 


Obs. Date 


MOS1 (ks) 


MOS2 (ks) 


PN (ks) 


0108260201 


July 01 2002 


62.44 


63.53 


43.68 


0110930701 


December 12 2002 


8.46 


10.05 


3.84 


0210270101 


December 19 2004 


25.42 


25.25 


23.79 


0210270201 


December 27 2004 


20.38 


19.92 


19.02 


0673310101* 


June 01 2011 


37.14 


42.96 


15.42 



2.1 Imaging Analysis 

Background-subtracted, exposure corrected images of all five 
pointings were produced in the soft energy band, 0.3 - 1.0 keV. 
The backgrounds for imaging an alysis were determin ed from cus- 
tomized blank-sky observations dCarter & Readll2007l) and closed- 
filter wheel observations to account for the instrumental back- 
ground. We use blank sky files with similar pointing directions and 
Galactic absorption column densities to M84 and M86. The closed- 
filter wheel counts images were subtracted from both the blank-sky 
and science exposures after renormalizing to match the total counts 
in each of these images outside of the field-of-view. The remaining 
blank sky emission was then subtracted from the science images. 
We also accounted for out-of-time (OOT) events in the PN detector 
images for each pointing using standard procedures. 

2.2 Thermodynamic Mapping 

Our XMM-Newton observations are sufficiently deep for high 
signal-to-noise spectra to be extracted from relatively small re- 
gions. This enables us to carry out detailed, spatially resolved 
measurements of the thermodynamic quantities across the en- 
tire mosaic. A ll spectral analysis was carried out using XSPEC 
dArnaudll2004l version 12.6). Instrumental background was deter- 
mined for the source spectra using closed filter wheel observa- 
tions of the detectors renormalized to have the same count rate 
in the 10 - 12keV energy band. We subtracted the out-of-time 
events from the PN spectra before performing these fits. We mod- 
eled the sky ba ckground for each spectral bin using the back- 
ground model of lUrban et al.l d201ll) . renormalized to account for 
the appropriate region geometries. This model assumes a sum of 
three compon ents: an absorbed power l aw from the unresolved 
point sources dDe Luca & Molendill2004l) . abso rbed thermal emis- 
sion w ith kT ~ 0.2 keV from the Galactic halo dKuntz & Snowdenl 
120001) and the unabsorbed thermal emission from the Local Hot 
Bubble dSidher et al]|l996l : iKuntz & Snowdenl l200Ch . Tests show 
that no additional background components are necessary. Allowing 
the temperatures, photon indexes, or normalizations of the back- 
ground components to vary did not change the results in a signifi- 
cant manner. 

2.2.1 Regions of Interest 

We identified the independent regions of interest for our spec- 
tral analysis by binning the background subtracted image into re- 
gions of equal signal to noise. Here we m ade use of the Weighte d 
Voronoi Tessellations binning algorithm bv lDiehl & Stalled d2006l) . 
which is a generalization of the Voronoi binning algorithm of 
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Figure 1. Background-subtracted, exposure corrected mosaic of five XMM-Newton pointings surrounding M86 (in units of ctss~'), utilizing the MOS+PN 
detectors in the energy band of 0.3 - 1 .0 keV. The four X-ray brightest Virgo Cluster members observed in this mosaic (M86, M84, NGC4388, and NGC 4438) 
are labeled. Diffuse gas is observed trailing to the northwest of M86, to the east of M86 in the direction of NGC 4438, to the south of M84, and to the east of 
NGC 4388. 



ICappellari & Copinl J2003h . We binned the image using two signal- 
to-noise thresholds (hereafter the low S/N and high S/N maps, re- 
spectively). Spectra were extracted from each spectral region, and 
appropriate response files were generated using the RMFGEN and 
ARFGEN tools. All spectral fits were performed in the energy band 
0.4-7. OkeV, although we ignored the energy range 1.2-1.8 keV in 
order to avoid the instrumental Al and Si emission lines in the MOS 
detectors. In the resulting energy band, there are roughly ~ 4000 



counts per region in our low S/N maps and 
gion for our high S/N maps. 



7000 counts per re- 



2.2.2 Spectral Modeling and Thermodynamic Quantities 

For the low S/N regions, we fit the spectra with an absorbed, sin- 
gle temperature AP EC plasma model. We utilized the photoele ctric 
absorption model of lBalucinska-Church & McCammoiil i ll 9921) and 
fixed the Galactic absorption column density to the value derived 
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i NGC 4388 
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Figure 2. An r-band image of the region surrounding M86 from the Sloan Digital Sky Survey. The X-ray surface brightness contours are overlaid in white, 
which show the stripped X-ray emitting gas trailing these galaxies. The stripped gas tail is especially pronounced for the M86 at the center of the image. Tails 
of stripped gas are observed to the south of M84 and to the northwest of M86. Diffuse soft emission is detected trailing eastward from M86 in the direction of 
NGC 4438. Faint, diffuse emission is also detected at the position of the edge-on spiral galaxy NGC 4402, located ~ 10' to the north of M86. 



bv lKalberla eTafl d2005h Nm = 2. 13 x 1 20 cm' 2 . We assume the so - 
lar element abundance ratios given bv lGrevesse & Sauvall dl998h . 
The temperature, metallicity, and the normalization were all inde- 
pendent free parameters. Best fit model parameters were obtained 
by minimizing the modified C-statistic in XSPEC. All uncertainties 
quoted are 68% (AC =1) confidence intervals. 

The high S/N spectra were fitted with a similar model, but 
with an additional thermal component to account for the ambient 
Virgo Cluster emission viewed in projection. The temperature of 



this thermal emission was fixed to 2.3 keV, in agreement with mea- 
surements of the di ffuse Virgo ICM at the same distance from M87 
dUrban et alj|201 lh . The metallicities of the two thermal compo- 
nents were linked to be the same. The exact choice of the tempera- 
ture of the ambient Virgo ICM does not have a strong influence on 
the measured temperature or metallicity for the second component. 

To search for multi-temperature structure in the high S/N map, 
we also attempted to fit the spectra with a two-temperature plus 
cooling flow model (APEC+APEC+MKCFLOW). For the cooling 
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flow component we set the upper and lower temperatures to 1 keV 
and 0. 1 ke V, respectively. All of our measured temperatures have 
an uncertainty of less than ~ 10%, regardless of the model. 



3 RESULTS 

The background subtracted, exposure corrected image of the region 
surrounding M86 is shown in Figure [TJ To the north-west of M86 
a long tai l of diffuse gas is observ ed, w hich has been discus sed in 
detail by iFinoguenov et alj d2004) and lRandall et al.l J2008h . The 
length of this tail is ~ 150 kpc in projection, although given the 
large line-of-sight velocity of M86 with respect to the ambient 
Virgo ICM this is likely a conservative lower limit to its true length. 
The width of this tail in projection, however, is roughly uniform 
along its entire length. To the east of M86, diffuse emission is ob- 
served out to the neighboring galaxy NGC 4438. Diffuse gas is also 
detected for a distance of approximately ~ 15 kpc to the south of 
M84. These tails of X-ray emission trailing M86 and M84 extend 
beyond the bulk of the optical galaxy emission, as shown in the 
Sloan Digital Sky Survey image of this field (Figure|2j. There is no 
evidence for diffuse emission that bridges between M84 and M86. 
There is also tentative evidence for diffuse X-ray emission trailing 
to the east of NGC 4388. 



3.1 Temperature Structure 

Maps of the temperature and metallicity structure across the re- 
gion are shown in Figure [3] Cool gas with a temperature of 
~ 1 keV is detected in various directions around M86, with the 
the majority of this gas phase located in the regions immedi- 
ately surrounding the galaxy and to the north-west of the galaxy. 
This ~ 1 keV component traces the morphology of M86's X-ray 
tail, and originates from the ram pressure stripping of M86 (e.g. 
Forman et al.lll979l: iFabian et al.lll980l: lNulsenlll982l : IWhite et all 
1 99 ll : IFinoguenov et al.ll2004lRandall et al.ll2008h 

To the east of M86 is a site of cooler X-ray emitting gas, ob- 
served to have a temperature of ~ 0.6 keV, spatially coincident 
with the brighter diffuse X-ray emission spanning between M86 
and NGC 4438. Gas with a temperature of ~ 0.6 keV is also ob- 
served just to the north of M86, extending to the northwest in 
the same direction as the ~ 1 keV component. Fitting each spec- 
tral bin in the region between M86 and NGC 4438 with a two- 
temperature plus cooling flow model results in no detection of sig- 
nificant cooling flows. For any particular bin, these fits constrain the 
spectroscopic cooling rate of X-ray gas from 1.0 keV to 0.1 keV to 

< 0.01M o yr~' , and the upper limit for the integrated spectroscopic 
cooling rate across all of the bins between these two galaxies is 

< 0.1M o yr». 

In the vicinity of M84, we detect the presence of cool (~ 
0.7 keV) gas trailing to the south of the galaxy, out to a distance of 
~ 15 kpc. This is significantly larger than the scale of AGN feed- 
back in this system (see Figure |4). The temperature of this gas is 
observed to increase smoothly to the south, until it becomes in- 
distinguishable from the ambient ICM. The presence of cool gas 
trailing to the south of M84 suggests ram pressure stripping of its 
X-ray halo amidst northward motion. 

We use the presence of this tail to provide an order-of- 
magnitude estimate of the rate at which X-ray gas is stripped from 
M84. In order to determine the gas mass in the tail, we assume it 
has a prolate ellipsoidal geometry with a semi-major axis of 7.5 kpc 



and a semi-minor axis of 4.2 kpc (shown in Figures|4j- This geom- 
etry is based on the extent to which we can detect cooler (~ 1 keV) 
gas to the south of the cavity regions. We only include the region 
to the south of the AGN-driven cavities seen in Chandra obser- 
vations of M84. Based on the surface brightness enhancement in 
this region, we estimate that the density of the gas in the tail is 
~ 0.005 cm -3 , roughly a facto r of 2 denser than the surrounding 
Virgo ICM dUrban et al.ll20l ]]). With these assumptions, we deter- 
mine that a mass of ~ 7 x 10 7 M o of X-ray gas has been stripped by 
the bulk motion of M84 through the Virgo ICM. Using the sound 
crossing time as an estimate for the time scale over which this strip- 
ping occurred and an assumed sound speed of 700 km s , the over- 
all mass stripping rate is calculated to be ~ 7M Q yr~' . Since there is 
no evidence for jumps in the ICM temperature or surface brightness 
indicative of a shock front to the north of M84, this galaxy is likely 
traversing through the ICM at sub-sonic speeds^]. Resultingly, the 
mass stripping rate we estimate is likely an upper limit as to the 
true rate at which M84's X-ray halo gas is being stripped by ram 
pressure. 

To the west of M86 in the direction of M84, a clear temper- 
ature gradient is observed over the entire region. Between the two 
galaxies is a ~ 20 kpc band of emission consistent with the ambient 
Virgo ICM, in both one- temperature and two- temperature fits. 



3.2 Metallicity Structure 

Although the single-temperature and two-temperature fits quali- 
tatively agree in their metallicity measurements (see Figures [3p 
and[3ji), we focus our discussion on the measurements from the 
two-temperature fits. These measurements are performed at higher 
signal-to-noise and are less sensitive to biases in metallicity mea- 
surements that aris e by modeli ng multi-phase emission with a 
single temperature l lBuote| [2002). For comparative purposes, the 
metallicity of the Virgo cluster at the same radius from M87 
dUrban et al.l201 ll) is approximately Z ~ 0.3 Z Q . 

The metal-richest gas is located near the center of M86, which 
has a measured metallicity of at least ~ 1Z Q . Metal-rich gas (> 
0.5Z o ) is observed in the direction of M86's stripped tail, out to dis- 
tances of ~ 90-100 kpc. There is no strong evidence for metal-rich 
gas spanning between M86 and NGC 4438. However, it is possible 
that a metallicity-measurement bias arising from unaccounted tem- 
perature phases is still present here. A smaller enrichment in metals 
(Z > 0.5Z o ) is observed near the core of M84. 



3.3 Ultraviolet Observations with GALEX 

We have also utilized publicly available ultraviolet (UV) observa- 
tions of the region between M86 and NGC 4438 with GALEX in or- 
der to search for star formation in this region. We used observation 
ID "NGA_Virgo_MOS10" (observation date 2004-03-11), part o f 
the GALEX Ultraviolet Virgo Cluster Survey dBoselli et al.lboi ll) . 
which is well centered between the two galaxies. There are no sig- 
inficant detections of young stellar populations in the far UV image 
of this field. Converting the observed s urface brightness of far-UV 
photons into a star formation rate using lKennicutj < fl998l) . we place 



1 Since M84 and M87 have similar line-of-sight velocities, it is also un- 
likely that M84 has a significant velocity component through the Virgo ICM 
along the line-of-sight. This maximizes the likelihood that such a shock 
front, if present, would be detected. 
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kT (Single Temp) 



Z (Single Temp) 




kT (Two Temp) 






0.302 0.837 0.936 0.989 1.04 1.1 1.17 1.23 1.28 1.4 1.7 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 



Z (Two Temp) 




0.125 0.825 0.905 0.942 1 1.06 1.08 1.16 1.31 2.02 5 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

Figure 3. Temperature and metallicity maps of the ~ 1° region surrounding M86, including the galaxies M84, NGC 4388, and NGC 4438. The upper row 
shows our low S/N (~ 4000 counts per bin) temperature and metallicity maps. The values show our best-fit parameters using a single-temperature thermal 
emission model. The lower row shows temperature and metallicity maps from our high S/N spectra (~ 7000 counts per bin). The values in these maps are 
derived from a two-temperature fit, with the second temperature fixed to 2.3 keV and the metallicities of both components tied to the same value. This two- 
temperature fit optimally identifies regions hosting gas significantly cooler than the ambient Virgo Cluster gas, which has an assumed temperature of 2.3 keV. 
The inclusion of an additional temperature component also reduces metallicity measurement biases. 



a strict upper limit on the star formation rate in the region of the Ho- 
filaments of < 0. 1M G yr~' . 



4 DISCUSSION 

Using a mosaic of XMM-Newton observations of the region sur- 
rounding M86, we have identified multiple sites of stripped X- 



ray emitting gas associated with all four galaxies: M86, M84, and 
NGC 4438, and NGC 4388. Although the evidence for stripped 
X-ray emitting gas near NGC 4388 is weak, radio observations 
of this galaxy show a tail of H I gas trailing to the north of 
this galaxy, with a velocity structure consistent with NGC 4388 
(e.g. lOosterloo & van Gorkoml l2005h . All four of these galaxies 
are therefore undergoing some form of ram pressure stripping, al- 
though each galaxy has its own unique observational signatures 
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associated with it. The stripped X-ray emitting gas trailing to the 
northwest of M86 and the H I tail trailing to the north of NGC 
4388 have been detected and discussed in previo us work (e.g. 
iRandall et al]|2008l ; lOosterloo & van Gorkomlbdoll) . The stripped 
X-ray emitting gas to the south of M84 and the cool X-ray emitting 
gas stretching from M86 to NGC 4438 have not been discussed at 
lengtrB We expand our discussion of these two sites of ram pres- 
sure stripping below. 

While all four of the galaxies lie near one another in projec- 
tion, the X-ray data combined with absolute velocity measurements 
along the line-of-sight suggest that only NGC 4438 and M86 are 
interacting with one another in any substantial way. There is no 
evidence in either the X-ray imaging or spectral maps for diffuse 
emission or cool gas bridging between any other pair of galaxies in 
a manner suggesting galaxy-galaxy interactio ns. The reported in- 
tracluster light (ICL) between M86 and M84 dRudick et ai]|201fj) 
therefore appears likely to be a projection effect. 



4.1 Correlation of the Coolest X-ray Gas and Ho- Filaments 

In Figure|5ji, we show the temperature map of Figure|3]; zoomed in 
to better show the ~ 0.6 keV gas plume to the east of M86. W e have 
overlaid the Ha filaments identified bv lKennev et"ai] I2008T) . which 
clearly show that the Ha filaments are located at the same positions 
as the lowest temperature phases of X-ray emitting gas. This region, 
at least th e half closer to M86, is also obs erved to be a site of HI 
emission (Oos terloo & van Gorkoml 12001 , although based on the 
velocity structure of that emission it may be primarily associated 
with NGC 4388 to the south of M86. On the other hand, some of 
the H I filaments detected ne ar M86 have velocity structure similar 
to M86 l lKennev et ai]|2008l) . 

A similar spatial coincidence between cool X-ray emitting gas 
and Ha filaments has been seen in a stripped tail behind the galax y 
ESO 137-001 falling into the cluster Abell 3627 JSun et al.ll2007h . 
The stripped gas trailing behind ESO 137-001 is a site of in situ 
star formation taking place outside of the galaxy within the ICM. 
A similar coincidence of cool X-ray emitting gas (at temperatures 
of ~ 0.6 keV) with filaments of Ha emission (at temperatures of 
~ 10 4 K) has also been observed around the BCGs of nearby cool 
core clusters such a s the P erseus Cluster dSanders & Fabian 20071; 



Fabian et al. l I2008L 1201 ll). Centaurus Cluster Zanders & Fabiarj 
2002Ll2008h , and M87 l lWerner et al.l2O10L |2012[) . 



4.2 On The Origin of the Ha Filaments 

The origin of this ~ 10 4 K phase of gas spatially coincident with 
soft X-ray emission between M86 and NGC 4438 is likely a previ- 
ous collision between the two. The existing multiwavelength data 
are all consistent with this scenario, although further data will be 
required to rigorously confirm this. The X-ray emitting gas in this 
filament region appears poorer in metals than the gas detected at 
the center of M86, implying that this cool X-ray emitting gas likely 
originated in NGC 4438. Measurement biases, due to the presence 
of multiple temperature phases unaccounted for in our models can- 
not be ruled out as the reason for the observed low metallicity in 
this region, however. 

In the absence of sources of heating, this stripped gas may 



be expected to cool to star-forming temper atures, as is obs erved 
in the galaxy ESO 137-001 in Abell 3627 JSun et alj|2007r) . Un- 
certainties in the filament geometry (specifically the depth of the 
filament regions along the line of sight) limit our ability to deter- 
mine the cooling time of the cool gas unambiguously. Based on 
the combined analysis of all of the spectral bins in this filament re- 
gion discussed above, we can place an upper limit on the integrated 
cooling rate from 1 keV down to 0. 1 keV across the whole region of 
< 0. 1M yr~' . This is in agreement with the upper limit on star for- 
mation imposed by far UV emission of this same region discussed 
above. 

Spiral galaxies such as NGC 4438 typically have X-ray halos 
with an observed temperature of ~ 0.2 keV, while the gas associ- 
ated with the group halo of M86 has a temperature of ~ 1 keV. 
The presence of ~ 0.6 - 0.7 keV gas in between the two galaxies 
suggests a possible mixture of gas phases. This "bridge" of cool 
gas between the two galaxies may arise either from ram pressure 
stripping of NGC 4438 by the surrounding Virgo ICM, or from 
the recent encounter between NGC 4438 and M86. The stripping 
medium in this latter scenario could either be the X-ray gas halo of 
M86 or the intragroup medium surrounding it, both of which could 
result in the observed filament of ~ 0.6 - 0.7 keV gas. 

The spatial coincidence between Ha emitting filaments with 
~ 0.6 keV X-ray emitting gas is similar to that observed in the cen- 
ters of cool core clusters undergoing AGN feedback. It is there- 
fore possible that similar physical processes drive the evolution of 
these distinct filame nt regions, at least after an initial disturbance . 
As discussed in (e.g. lFabian e"tai]|2008ll2oTll ; IWerner et al.ll2012h . 
a potential energy source powering the FIR through UV emission 
(including the Ha+N II optical line emission) in cooling core clus- 
ters are hot ICM particles penetrating into the filamen ts of cold gas, 
result i ng in collisional ion i zation and heating (see iFerland et al.l 
120081 l2009h - 1 Wern er et al. I d2012l) propose that shear instabilities 
due to shocks or the motion of the filaments through the ambient 
ICM may facilitate mixing between the hot and cold gas phases. 
The physics of the multiphase gas in the M86 filament region, 
which is exposed to shear instabilities as it moves through the Virgo 
ICM and merger shocks due to the infall of the M86 group, may 
therefore share similarities with the physics of the optical emission 
line nebulae seen in the cooling cores of galaxy clusters. One key 
difference between the filament region discussed here with those 
observed near cool core clusters, however, is the nature of the initial 
disturbance. Rather than instabilities triggered by AGN feedback, 
the collision between M86 and NGC 4438 would be the source of 
the shearing here. Such triggering may be further encouraged by 
M86's supersonic motion and subsequent shock heating. 

Sites within ~ lOkpc of the two galaxies have recently been 
shown to host CO(l — » 0) line emission, suggesting that molecular 
hydrogen gas is presen t between the two galaxies for time scales 
of at least ~ 100 Myr jPasvra et al.ll2012T) . This molecular gas is 
inferred to originate primarily from NGC 4438, similar to what we 
infer for the X-ray emission. Further observations that place stricter 
constraints on the star formation rate and detailed spectral proper- 
ties of this region, as well as the mass of cold gas present from 
emission lines such as C II, will be required to test the connection 
between the Ha emission presented here and that observed in cool 
core clusters. 



2 The diffuse tail of gas trailing to the south of M84 was discussed briefly in 
IRandall et al. 1 2008), although there were insufficient counts in the Chandra 
data to resolve its thermodynamic structure. 



4.3 Combining AGN Feedback with Ram Pressure Stripping 

The detection of a clear tail of X-ray emission trailing to the south 
of M84 offers another example of a galaxy where feedback and ram 
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Figure 4. A zoomed in view of M84 with XMM-Newton. Both of these images show the same field of view. The dashed ellipse denotes the region which 
we use to estimate the rate at which X-ray halo gas is being stripped from M84. Left: The background subtracted, exposure corrected image of M84 in the 
0.3 - 1.0 keV energy ba nd. The black contours are derived from Chandra observations of this galaxy, which show the extent of the AGN-driven cavities 
iFinoguenov et alj2008i) . Right: The corresponding temperature map of this same region, in units of keV. The bin map is similar to that described in $2.2. II 
although with a only ~ 250 net counts per bin in order to resolve smaller scale features. Each bin is fit with a single-temperature plasma model with Galactic 
absorption (PHABS X APEC). The same Chandra contours are overlaid in black. The brighter gas trailing to the south of M84 is shown to be significantly 
cooler than the ambient ICM. 



pressure stripping work in conjunction to strip its X-ray halo. Ram 
pressure strippin g combined with AG N feedback are also observed 
in Sersic 159 -03 ( Werner et al.l201 1| ) and the Virgo elliptical galaxy 
NGC 4552 dMachacek et alj|2006j|bl) . and similar processes are 
obser ved at larger scales in the Ophiuchus Cluster dMillion et al .1 
l20ld) and MACS J 193 1.8-2634 dEhlert et al.|[201 lh . In these sys- 
tems, the combined influences of these two processes can rapidly 
suppress or even destroy cool cores associated with the central re- 
gions of the galaxy/galaxy cluster. It appears as though a configu- 
ration where both AGN feedback and ram pressure stripping occur 
simultaneously is more common than initially expected. 

The energy injected into the X-ray halo from AGN feedback 
may allow M84 to be stripped more efficiently now than in its past. 
At its current stripping rate of < 7M G yr~' estimated above, all of 
the X-ray gas should be stripped from this galaxy within ~ 200 - 
300 Myr, a time scale likely far shorter than the time elapsed since 
M84 first crossed the Virgo Cluster's virial radius. We determine 
whether losses from stellar winds may be sufficient to replenish 
the X-ray halo of M84, assuming a single-age passively evolving 
stellar po pulation wit h a Sal peter IMF. We use the stellar mass-loss 
rate from lCiotti et al. 

icon 

to estimate this rate, M+(t), as 



M„(t) - 1.5 x IO^'Le^ 13 



(1) 



where Lb is the present-day B-band luminosity in units of Lb.o 
(L b = 4.5 x 10 10 L B , o for M84; iFinoguenov & Jonesll2002l) and h s 
is the age of the stellar population in units of 15 Gy3~The local 
stellar mass-loss rate, assuming that the current stellar population 



This formula is valid in the range from 0.5 to over 15 Gyr. 



is 10 Gyr old, is of the same order-of-magnitude as the rate of ram 
pressure stripping (M 4 ~ lM yr~'). 



4.4 On the Frequency of Stripping in Cluster Member 
Galaxies 

Perhaps surprisingly, each of the four X-ray emitting galaxies de- 
tected in this mosaic show evidence of being stripped as they tra- 
verse through the ICM. Although these galaxies share some com- 
mon features, no two galaxies are observed to have stripped gas 
tails sharing all of the same properties. The prevalence of gas trail- 
ing in the wake of galaxies in this busy region of the Virgo Cluster 
suggests that stripping processes may play an important role for 
nearly all galaxies in clusters. The particular observational signa- 
tures associated with the stripping process may depend strongly on 
the morphology, mass, and gas content of the host galaxy, but it 
appears typical that gas initially hosted by galaxies is efficiently 
stripped by the ICM. 

For more massive elliptical galaxies hosting large masses of 
X-ray emitting gas like M86 (likely the central galaxy of an in- 
falling group) and M84, the stripping is most readily observed as 
tails of X-ray emitting gas trailing the galaxies. The stripped X-ray 
gas behind these two galaxies, however, suggest that the astrophys- 
ical processes that drive the stripping may vary between different 
sites. One key difference between the stripped X-ray gas between 
M86 and M84 is the relative length of the apparent tails, which may 
indicate differences in the gas physics between the two sites. The 
tail behind M86 remains roughly uniform in width and coherent 
for distances of at least ~ 150 kpc in projection without significant 
evidence for mixing with the surrounding ICM. The true length of 
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this tail is also likely significantly larger than its observed length 
in projection due to M86's high velocity with respect to Virgo 
(~ 1500 km s -1 ) along the line-of-sight. Maintaining such cohesion 
over these distances suggests that the stripped gas is undergoing a 
laminar flow even as the galaxy traverses through the ICM super- 
sonically, and subsequently suggests that the stripped gas may have 
high viscosity. The gas trailing behind M84, on the other hand, is 
only observable for distances of ~ 15 kpc in surface brightness and 
temperature maps, and the lower velocity of M84 with respect to 
the Virgo ICM along the line-of-sight (~ 200 - 300 km s -1 ) indi- 
cates that this tail is probably not significantly longer than its pro- 
jected length. The temperature structure in this tail is also observed 
to increase smoothly with distance from M84. These may be evi- 
dence for more efficient mixing between the stripped gas and the 
surrounding ICM at this site, perhaps associated with a more turbu- 
lent flow of the stripped gas. Deeper observations will be required 
to determine the full length of the stripped gas trailing M84 and 
subsequently the extent to which turbulence and magnetic fields 
may be necessary to account for the thermodynamic structure in 
M84's tail. 

Spiral galaxies contain significant masses of gas at much 
lower temperatures (~ 10 - 100 K). The observed differences in the 
stripped cold gas behind NGC 4388 and NGC 4438 suggest that the 
manner in which these phases interact and intersperse with the ICM 
may vary also significantly from site to site, and that the presence 
of shear instabilities or shocks may play an important role in these 
processes. This also appears to be true for the X-ray halos, espe- 
cially when the disruptions due to the surrounding environment are 
more extreme. Deeper multiwavelength observations will again be 
necessary to fully understand the range of gas temperatures present 
in these stripped tails and test the stripping processes from start to 
finish. 



5 CONCLUSIONS 

We have measured the X-ray spectral properties in a contiguous 
~ 1° x 1° region spanning four Virgo Cluster galaxies, allowing us 
to view interactions between them. Our key results are as follows: 

1) There is clear evidence for cool (~ 0.6 keV) gas trailing to 
the northwest of M86 as reported in previous studies, to the east 
of M86 in the direction of NGC 4438, and to the south of M84. 
M86 appears to be surrounded in a halo of ~ 1 keV gas, which we 
interpret as the intragroup medium of a recently merged subc luster 
undergoing ram pressure stripping by the ambient Virgo ICM. We 
also observe tentative evidence for stripped X-ray gas to the east of 
NGC 4388. 

2) The presence of stripped gas associated with all four of the 
largest galaxies in this field of view indicates that, while the partic- 
ular astrophysical details may vary from site to site, ram pressure 
stripping is a common occurrence in cluster member galaxies, and 
has important implications for the heating and enrichment of the 
ICM. 

3) Filaments of Ho- emission surrounding M86 are observed 
at locations spatially coincident with enhanced X-ray emission at 
temperatures of ~ 0.6 keV. These filament regions have similar X- 
ray spectra to Ho- filaments observed in the centers of cool core 
clusters. One key difference between this filament region and those 
detected in cool core clusters, however, is the likely source of the 
initial disturbance that triggers their formation. 

4) AGN feedback and ram pressure are both observed in M84, 
and in tandem these two processes may be encouraging particu- 



larly efficient stripping. Without replenishment, M84's X-ray halo 
gas would be depleted on time scales of ~ 300 Myr. Mass loss 
from stellar winds, however, appears to be sufficient to replenish 
the M84's X-ray halo, at least in principle. 

5) Our results for M86 and M84 suggest that there is no signif- 
icant interactions taking place between them. Of these four galax- 
ies, only M86 and NGC 4438 appear to be undergoing any notice- 
able galaxy-galaxy interactions. 

Disentangling the interactions of the infalling M86 galaxy 
group in the directions of M84, NGC 4338, and NGC 4438 will 
require further observations at multiple wavelengths. 
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Figure 5. The correlation between the cool X-ray gas and Ha emission detected by Kenney et al. 2008. (a) Temperature map of Figure [3Ji zoomed in to 
emphasize the temperature structure surrounding M86 and overlaid with contours of Ha imaging. There is a clear correlation of the Ho- emission with the 
coolest gas phases in both the eastern and northern plumes extending out from M86. The centers of M86 and NGC 4438 are denoted by the cyan crosses, (b) 
The X-ray surface brightness image of Figure[T] zoomed in to the same region as (a). The Ha image contours are overlaid in cyan. 
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